Analysis of constrained theories without use of 
primary constraints. 
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Abstract 

' It is shown that the Dirac approach to Hamiltonization of singular the- 

ories can be slightly modified in such a way that primary Dirac constraints 
£""•«» ■ do not appear in the process. According to the modified scheme, Hamilto- 

OO \ nian formulation of singular theory is first order Lagrangian formulation, 

further rewritten in special coordinates. 

O 

in 

1 Introduction 

Let us consider Lagrangian theory with the action and the equations 
of motion being 



X 

c3 ■ fdLYdL 



S = JdrL(q A ,q A ), A = 1, 2, • • • [A], (1) 

0. (2) 



\dq A J dq A 
It is supposed that the theory is singular 

rank ww = [ ' ] < [A] ' (3) 

According to Dirac [1], Hamiltonian formulation of the theory is 
obtained as follow. First step of Hamiltonization procedure is to 
define equations for the momenta pa- Pa = Jpr • The latter can be 
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considered as algebraic equations for determining of velocities q A . 
According to the rank condition (jSJ), [i] equations can be resolved for 
q l and then substituted into the remaining ones. By construction, 
the resulting equations do not depend on q A and are called primary 
constraints $> a (q,p) of the Hamiltonian formulation. The equations 
Pa = Jpr are then equivalent to the following system 

q l = v\q A , Pi ,q a ), (4) 

<$>a(q A ,Pa,Pi) = 0. (5) 

On second step one introduces an extended phase space with the 
coordinates (q A ,PA, v a)- By definition, Hamiltonian formulation of 
the theory (0) is the following system of equations on this space 

q A = {q A ,H}, p A = {p A ,H}, 

$a(g A ,p B ) = 0, (6) 

where { , } is the Poisson bracket, and it was denoted 

H(q A ,p A ,v a ) = Ho(q A , Pj )+v a ^ a (q A ,p B ), (7) 



H = p^-L + q ' 



dq c 



(8) 

q 1 ^v 1 (q A ,Pj,q a ) 



It is known [2] that the formulations (J2J and © are equivalent. 

Aim of this part of the Dirac procedure is, in fact, to represent 
system of differential equations (J2J) in the normal form (with deriva- 
tives of variables separated on the l.h.s. of equations) . Adopting this 
point of view, we demonstrate below that Eq.(j5J), as well as the vari- 
ables not necessary in the process. In Sect. 2 we formulate 
Hamiltonization procedure for the theory without use of the 
momenta p a . As a consequence, the primary constraints (0) do not 
appear in the formulation. Concerning to Eq.p|). it has simple in- 
terpretation [3] as (invertible) change of variables on configuration- 
velocity space 1 . Expression for the change will be obtained here 
as sufficient condition for conversion of equations of motion (see 
Eq. (jl7|) below) into the normal form. Equivalence with the Dirac 
procedure is shown in Sect. 3. Some useful identities are collected in 
the Appendix. 

1 In this respect, the formulation presented below is more close to the case of nonsingular 
theory, where transition from the Lagrangian to the Hamiltonian formalism is change of 
variables. 
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2 Hamiltonization without use of primary con- 
straints 



We present here modification of standard Hamiltonization procedure 
[1,2] which do not implies appearance of the primary constraints in 
the formalism. 

1 Step. System of differential equations of second order (J2J) can 
be rewritten in a first order form. To achieve this, one introduces 
configuration-velocity space with the coordinates (q A ,v A ). The ini- 
tial system and the following one are equivalent 2 



q A = v A , M AB v B = K A . (9) 
Here it was denoted 

dv A dv B ' ' dq A dv A dq B 



/ \ d 2 L(q,v) - . . dL d L R 
MAB(q,v) = —^-, K A (q,v)^—j-—— w v B . (10) 



We use symbols with bar to denote functions on the configuration- 
velocity space: A = A(q A ,v B ), and symbols without bar for func- 
tions on extended phase space (see below) 

A = A(q A , Pj ,v a ) = A(q A ,v a ,v% i(q A p ., va) . (11) 

Without loss of generality [2], we suppose that the variables q A has 
been enumerated in a way adjusted with Eq.flSJ), then 

MAB= \M a] M a /3 j 
A=(i,a), i = l,2,-.-[i], a=l,2,..-[a] = [i4]-[i], 

rankM AB = [i], detM^ 0. (12) 

- —ij — — jk , 

Inverse matrix for My is denoted M , one has .\ /, ; .\ I = Si . 
2 Step. According to the rank condition, [i] equations among 

(J9J) can be resolved in relation of if: v l = M (Kj — Mj a v a ) and 
then substituted into remaining [a] equations with the result being 

[M a/3 - M ai M" 'Mjpjv 13 = K a - M ai M 3 Kj. It must be rank[M af3 - 

- — ij - 

M ai M Mjp] = (if not, one would be able to resolve the equations 

2 The systems J5J and J5J arc equivalent in the usual sense [2]: a) if <7o( T )i ^o( T ) i s solution 
of then qo(i~) is solution of J3Ji b) if 90 (t) is solution of then qo(r), po(r) = qo is 
solution of 151 . 
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in relation to some of v a , in contradiction with @). Then one has 
the identity 

M a(i - M ai M 3 Mjp = 0. (13) 
As a consequence, the following quantities 

C a A = (-M a fi ] \5jy (14) 
form a basis on space of null-vectors of the matrix M A b (see also 



C a A M AB = 0. (15) 
In the result, our equations © acquire the equivalent form 

q A = v A , (16) 



if = M 3 (Kj — M ja i) a , ) (17) 



K a - M ai M Kj = C a A K A = 0, (18) 

where the system ()18|) is an algebraic one, (also, it may contain 
functionally dependent equations). 

3 Step. Eq.(|17J) is not in the normal form. Then one asks 
on existence of special coordinates (q A ,v A ) on (q A , v A )-space which 
reduce the system to the normal form. We demonstrate here the 
following statement. 

For singular theory (|T|)-(|3|) there is exist invertible change of vari- 
ables of the form 

a \ 

v i (q A ,p j ,v a ) 7 rank— = [^, (19) 
d Pj 

which reduces the system (|17j) to the normal form. 

To make a proof, let us write Eq. (jl7|) in terms of new variables 
(see our notations 

dv 1 ( BtP \ difl a 

MtJ dp~ k Pk + { Mij d^ + Mia ) ia = Ki ~ m 



1 


q A \ 


fq A \ 












V 


v i J 


{ Pi J 





4 



It reduces to the normal form if the functions v 1 obey the equations 
(note that My is the only invertible matrix at hand) 



dv 1 
d Pj 



M l K 



dv 1 



-M' J M 



ja- 



(21) 



These equations determine v 1 modulo of an arbitrary function A(q A ), 
the latter will be omitted in what follow. Suppose that there is exist 
a solution of Eg .([2 1)1 . and let pi(q A ,v A ) be inverse transformation. 
Expression for the latter can be easily obtained. Actually, using 
Eq.(|2*T|l. derivatives of the identity v z (q A ,Pi(q J 



v 1 



be written in the form 



dpi 



M*. 



d 2 L 



M, 



= v can 

d 2 L 



Solution of the system gives expression for v l in implicit form 



Pi 



dL(q A ,v\v a ) 
dv 1 



(22) 



dL 



It is necessary condition which turns out to be a sufficient condition 
also. Actually, Eq.(|2*2*|) is invertible due to Eq.®, and defines some 
function v t (q A ,Pj,v a ). Then derivatives of the identity Pi 
imply immediately the equations (|21j) (see also Appendix). 

Curiously enough, the reasoning works for non singular theory 
also. To find normal form for non singular equations MabV B = Ka, 
one applies the inverse matrix M ab- Equivalently, one can search 
for a change v A (q A ,pA) which reduces the system to the normal 
form. It gives the standard Hamiltonian formulation of the theory. 

4 Step. Let us rewrite the equations (fT(3|) -(fT8 |) in the coordi- 
nates (fTTJj) . Properties of the transformation v l = v l (q A ,pj, v a ) are 
presented in the Appendix. In particular, it turns out to be linear 
on v a , and has the representation 



v i (q A ,p j ,v a ) 



OH, 



dpi 



where 



Hr = Ho(q A , pj) — f a (q A , Pj)i 



fa 



dL 



(23) 



(24) 



and Hq coincides with the Hamiltonian (JHJ) of the Dirac scheme. 
Moreover, in the coordinates all the quantities of the first order 
Lagrangian formulation can be rewritten in terms of the functions 
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H and f a , see Appendix. Using Ens.(fT5|). l|23|» . (I2H). pE)) . (EE2l. the 
system ()16p -(|18) l can be identically rewritten in the form 

q a = v a , (25) 

q = W Pl = ~w (26) 

A Q/3 (q A ,Pjy + H a (q A , Pj ) = 0. (27) 
Here it was denoted 

A a/3 = {/ Q , fp} R - (-^ - , H a = -{f a , H }r - -q^>( 28 ) 

and { , }r is reduced Poisson bracket (i.e. the Poisson bracket on 
q\Pi- subspace). For any functions A, B of the variables (q A ,pi, v a ) 
it is defined by 

r . D1 dAdB dAdB 

It finishes first stage of the Hamiltonization procedure. Differential 
equations are presented in the normal form 1)26)1 . One notes also 
that algebraic equations ()27j) turn out to be linear on v a . They co- 
incide with the Dirac second stage system (for revealing of secondary 
constraints): H} = (see also the next section). 

In resume, initial equations of motion (0) for singular theory 
have been rewritten in special coordinates ()19)1 on the configuration- 
velocity space. Motivation for making of the change was to find 
normal form of differential equations ([17)1 . It turns out to be suffi- 
cient to transform the variables v % only, defining the new variables 
Piiyi) in accordance with Eq.()22|). One concludes that representa- 
tion of singular Lagrangian dynamics in the normal form ()25j) -(l27 |) 
do not implies any extension of the configuration-velocity space. In 
particular, the primary Dirac constraints do not appear in the for- 
mulation. We have demonstrated also that the Dirac second stage 
system {$ a , H} = represents part of Lagrangian equations of mo- 
tion (J18)l . rewritten in the coordinates (q A ,pi,v a ). By construction, 
the systems (0) and (j2~5|) - (j2~7j) are equivalent. It would be interesting 
to analyze local symmetries [4, 3] in this framework, as well as to 
compare the reduced scheme with other modifications of the Dirac 
procedure [5, 6]. 
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3 Equivalence with the Dirac procedure 



We discuss here relation among the system (j25|) - (|27jl and the stan- 
dard Hamiltonian formulation (jUJ). Starting from (J25j) - (|27|) . let us 
define an equivalent system on extended space with the coordinates 
(q , pj,v a ,p a ). By definition, (q A ,Pj, f a )-sector obeys the Eqs.((2SJ)- 
(|27|). while for p a one writes the primary Dirac constraints 

$a =Pa -fa(q A , Pj) = 0. (30) 

Then Eq.([27|) can be rewritten identically in the form H} = 0, 
while the following equation: 

dH . , 

P- = ( 31 ) 

is a consequence of Eqs. (j25|) - (j2"Tj) . (JHOJ), and thus can be added to this 
system. After that, Eq. (j2"7j) turns out to be consequence of others, 
and can be omitted. The resulting system (|2*5j) . (|26|). (|3(Jj) . (j3T|) 
coincides exactly with the Hamiltonian formulation (jHJ) (modulo of 
notations). 

This discussion reveals that the only role, playing by p a in the 
Dirac formulation, is to represent equations of motion in completely 
symmetric form, with the Poisson bracket defined in relation to all 
variables q A ,PA- The momenta p a are, in fact, "Hamiltonian ghosts" 
of singular theory. The symmetrization seems to be unnecessary 
ingredient of this stage. Further stages of the Dirac procedure (on 
revealing of higher-stage constraints) do not require the symmetric 
form of equations. It may be more convenient to make further steps 
by using of the reduced formulation, and then to symmetrize the 
resulting equations. 



4 Conclusion 

In this work it was demonstrated that Hamiltonian formulation 
for singular theory (JT |) -lj3* j) can be considered as the first order La- 
grangian formulation Q on the space (q , v\ v a ), further rewritten 
in special coordinates (q A ,pi,v a ) of this space. The corresponding 
change of variables was find from the requirement that equations 
(jHJ) acquire normal form in the new coordinate system. The trans- 
formation v t (q A , pj , v a ) , defined implicitly by Eq. (J22|) . turns out to 
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be sufficient condition to this end. Search for the normal form do 
not implies use of an additional variables, as well as the primary 
Dirac constraints. The resulting formulation ()25j) - (j27j) involves al- 
gebraic equations (|27|) linear on v a . They coincide with the second 
stage equations of the Dirac procedure: {§ a ,H} = 0. The system 
(j2*B]) -(P7 j) and the standard one (JUJ) are equivalent. 
For the scheme presented one writes, schematically 

q A ^(q A ,v\v a )~(q A ,p u v a ). (32) 

It can be compared with the standard approach: according to [2], 
Hamiltonian formalism involves more variables 

q A ^ (q A ,v\v a ) -> (q A ,v\v a , Pl ,p a ). (33) 

Dynamics of (q A , w j4 )-sector is governed by the Lagrangian first order 
equations Q, while pa = -§^k- 
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6 Appendix 

We present here properties of the function v l (q A ,pj, v a ) defined im- 
plicitly by Eq. (}2"2"|) . as well as properties of the Lagrangian L in 
terms of the coordinates (fTTJj) : L(q A ,pj,v a ) = L(q A , v\ v a )\ v i . In 
particular, both functions turn out to be at most linear on v a in any 
singular theory (JD)-©. 

Derivatives of the identity pi = dL<yq d ^ - - \ v i {q A , Pj ,v a ) gi ye the re- 
lations (see our notations (JTTJ) , (JHJ)) 



dpj dv a Jon dq A dvidq A 



(34) 



Then the identity (|T3j) acquires the form (^lu*) — 0, so the 
quantity 



U(q A ,Pj) = 



(35) 
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does not depend on v a . One finds d ^ a<y Q p ,p ^ = M^Mj a , so the 

quantity M lJ Mj a does not depend on v a also. Moreover, Eq. (|33J) 
can be identically rewritten in terms of L 

fa(9 A i Pj) = ~(piV i -L) i (36) 



dv 

which implies that the quantity 

H R (q A , Pj ,v a ) = (p i v i -L), (37) 

is at most linear on v a . Integral of Eq. ()36|) gives the following ex- 
pression for H R 

H R (q A , Pj ,v a ) = H (q A , Pj )-v a f a (q A , Pj ), (38) 

where Hq does not depend on v a . Expression for this quantity 
in terms of initial Lagrangian can be obtained by comparison of 
Eqs.flHZj), (HP 



(39) 



Note that it coincides exactly with the Hamiltonian H (jSJ) of the 
standard scheme. From Eq. (}3Tj) one obtains useful relations 



dH R = vi dH R _ _dl_ 
d Pi dq A dq A 



flTT 

-^ = -U<l A ,Pi)- (40) 



In particular, the first equation states that our function v % is at most 
linear on v a and has the representation 

it A a\ ^ fa a 

v(q ,„„„) = _-_„. (41) 



It implies (see Eq. f!H4j) ) that the matrix M* J is at most linear on v a , 
and has the representation M ij = jj^. At last, Eqs.(J37J), (IHH|) . 
(}4~T|) imply that L is at most linear on v a and has the representation 

L(q A , Pj ,v a )= Pl ^-H R . (42) 

In resume, we have demonstrated, in fact, that in the coordinates 
(q a , P i,v a ) all quantities of the first order Lagrangian formulation 
(P)- fTUj) can be rewritten in terms of the functions H (q A ,pj), f a (q A , P j)- 
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